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Abstract. Equations of state for strange quark matter in beta equilibrium at high densities 
are used to investigate the structure (mass and radius) of compact objects. The chromodi- 
electric model is used as a general framework for the quark interactions, which are mediated 
by chiral mesons, a and if, and by a confining chiral singlet dynamical field, x- Using a quar- 
tic potential for x, two equations of state for the same set of model parameters are obtained, 
one with a minimum at around the nuclear matter density po and the other at p~5po. 
Using the latter equation of state in the Tolman-Oppcnhcimcr-Volkoff equations we found 
solutions corresponding to compact objects with R ~ 5 — 8 km and M ~ Mq. The phe- 
nomenology of recently discovered X-ray sources is compatible with the type of quark stars 
that we have obtained. 



INTRODUCTION 

Various effective models using quarks as fundamental dynamical fields, originally 
designed for the nucleon, have also been successfully used to describe infinite 
quark matter, and the resulting equations of state (EOS's) have been applied to 
investigate the structure of compact stars [1-6]. 

The chromodielectric model (CDM) [7-9], for example, provides a reasonable 
phenomenology for the nucleon [10, 11] and also allows us to obtain EOS's for 
dense quark matter. The model yields soliton solutions representing single baryons 
with three quarks dynamically confined by a scalar field, x> whose quanta can 
be assigned to ++ glueballs. When it is applied to quark matter in two or three 
flavors [12-14] the resulting EOS turns out to be relatively soft at large densities. 

Using a quadratic potential for the x field, Drago et al. [15, 16] applied the CDM 
to describe the inner part of neutron stars, obtaining masses in the range 1 — 2M 
and radii of the order 10 km or higher, with a small hadron crust of 2 km. In this 
work we consider an extension of the model used in Ref. [15], taking quartic instead 
of quadratic potentials. In addition to the structures found by Drago et al., the 
quartic model predicts another type of compact objects made out of quarks only, 
smaller and denser than neutron stars. 
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From the observational point of view, the recent discovery of X-ray sources, by 
the Hubble and Chandra telescopes, increased the plausibility that these sources 
might be strange quark stars [17-19]. In particular, the isolated compact object 
RX J1856. 5-3574 with a small radius does not show evidence of spectral lines or 
edge features [20, 21], reinforcing the conjecture for the existence of stars made out 
of strange matter. The phenomenology of these objects seems to be compatible 
with the small and dense quark stars reported in this work. 



THE MODEL 
CDM Lagrangian 

We write the CDM Lagrangian in the form [8-10] 

£ = £+£+£ +£ (I) 

where 

C q = i^d^ , = Idfi&a + \djt ■ 8^ - W{%, a) , (2) 

and W(n, a) is the Mexican hat potential. In the u, d sector the quark-meson 
interaction is described by 

£ g -meson = + IT ■ 7T-f 5 )ljj . (3) 

The last term in (1) contains the kinetic and the potential piece for the field: 

c x = ±d li xd»x-u(x). (4) 

The potential term is 



u(x) = \m\ 2 




6?7 4 \ x 2 



7 2 J (7M) 2 



(5) 



where M is the x mass. The parameterization used in (5) allows for a physically 
meaningful interpretation of the parameters 7 and r]: U (%) has a global minimum 
at x — and a local one at x — jM, and U{^M) = r) 4 M 4 (see Fig. 1). The height 
of the local minimum, B = (r/M) 4 ", is interpreted as a "bag pressure" and this is 
used to fix the parameters in U(x)- Assuming the wide range 0.150 <B 1 / 4 < 0.250 
GeV, one has 0.08 <i]< 0.15, using M = 1.7 GeV. We note that 7 is not a free 
parameter since the quartic term of U(x) must be positive and the cubic term 
negative, which implies 7 2 > 6^ 4 . In the soliton sector of the model, best nucleon 
properties are obtained for G = \J gM ~ 0.2 GeV (only G matters for the nucleon 
properties) and we keep that G in our quark matter calculations. 

In order to study strange quark matter, we add to the interaction Lagrangian 
(3) the term [22] 

^s— meson Ipslps ■ (6) 

X 

accounting for the coupling between the strange quark and the x held. 
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FIGURE 1. Quartic potential, Eq. (5), for fixed M and 7, in dependence of For comparison, 
the quadratic potential (7 — > 00) is also shown. The value of the potential at the second minimum 
can be interpreted as a bag pressure, B. 



Strange quark matter 

In order to study strange quark matter in beta equilibrium, an electron gas must 
also be considered. The mean-field energy per unit volume for strange quark matter 
in the CDM (plus electrons) is given by 




where the first two terms refer to quarks and the third one to the electrons, all 
described by plane waves. The degeneracy factor is a = 6 (for spin and color). The 
last term is the Mexican hat potential (with 7? = and f n = 93 MeV). The ki in 
(7) are the Fermi momenta of quarks and electrons. 

The quark masses in (7) are [22]: m u ^ = g u ,d& / {xfir) an d Tn s = 9s/x with the 
coupling constants given by g u = g (/ w +f 3 ), g d = g (f w - f 3 ) and g s = g (2f k - / w ) 
[f 3 = -0.75 MeV, f K = 113 MeV]. 

A variational principle applied to the energy density, Eq. (7), leads to two gap 
equations for a and X- 111 the interior of a compact star the matter should satisfy 
both the electrical charge neutrality and chemical equilibrium. These conditions 
should supplement the gap equations, and altogether we have a system of six 



Small quark stars in the chromodielectric model 



February 1, 2008 



3 



algebraic equations to solve at each baryon density p = (p u + Pd + Ps) /3 [here, pi = 
akf/(6n) 2 stand for each flavor density]. The solution of the system of equations 
are the meson fields, a and Xi an d the Fermi momenta, k u , kd, k s and k e . For the 
same set of model parameters we found two stable solutions, hereafter denoted by 
I and II, shown in Fig. 2. 
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FIGURE 2. Fermi momenta (solid lines) and quark masses (dashed lines) for solution I and 
II in dependence of the baryon density. For solution II the quark masses and the electron Fermi 
momentum almost vanish. The model parameters are g — 0.023 GeV, M — 1.7 GcV, 7 = 0.2 and 
rj = 0.12. The vertical scales are in fm _1 . 

For both solutions o is always close to f n . In solution I, the x fi e ld is a slowly 
increasing function of the density, remaining always smaller than ~ 0.05 GeV. For 
such a small x, the quartic potential and the quadratic potential are indistinguish- 
able (see Fig. 1), thus, in practice, solution I corresponds to the one obtained and 
used by Drago et al. [15] in the framework of the quadratic potential. Due to the 
smallness of the x field, quark masses are large and the system is in a chiral broken 
phase. The solution II exhibits a large confining field, x ~ iM (local minimum 
of U), independent of the density. The resulting quark masses are similar for the 
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three flavors and very close to zero (chiral restored phase). Therefore, the chemi- 
cal potentials in solution II are dominated by the Fermi momentum contribution, 
ji u ~ p d ~ p s and p e ~ 0, i.e. in solution II there are almost no electrons. Besides so- 
lutions I and II, there is an additional unstable solution corresponding to x ~ •yM/2 
[local maximum of U(x)]- 



Equations of state 

The energy per baryon number as a function of the baryon density (EOS) is 
readily evaluated for each solution. We stress that EOS-I is not sensitive to 7 and 
77 (since x is small), just depends on G. The saturation density occurs at a relatively 
low density (p = 1.2 — 1.6 po) and the behavior of EOS-I at intermediate densities, is 
similar to the hadronic EOS's (see Ref. [13] for the two flavors sector). The EOS-II 
is also insensitive to 7, but does depend on r\ [in fact, the dependence is on (r^M) 4 , 
as we have already pointed out]: the energy per baryon number increases with rj 
and so does the saturation density. For 77 ~ 0.12 the saturation density is p ~ 5po 
and the energy per baryon number is some 20% higher than for solution I at its 
saturation density. 




e / GeV fm" 3 e / QeV fm"' 

FIGURE 3. Pressure versus energy density for the two types of solution for various parameters 
77. ] The other model parameters are M = 1.7 GeV, g = 0.023 GeV and 7 = 0.2. 
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In Fig. 3 we present the EOS's in the form of pressure P = p 2 f^(j) versus 
energy density plots for solutions I and II and various rj. It is worth noticing that 
our results with two distinct EOS predicted by the CDM, are consistent with 
the results from perturbative QCD: it is remarkable indeed that, for 77 ~ 0.12 the 
CDM reproduces accurately the EOS's recently obtained in a perturbative QCD 
calculation (compare our third panel of Fig. 3 with figure 1 of Ref. [23]). 

Regarding energetics, both phases are almost degenerated at high densities and 
have similar shapes in the Pxe plane even at intermediate densities (or energy 
densities) in the narrow range 0.1 < rj < 0.12. In that region of 7], one solution is 
not clearly lower in energy than the other. However, we should point out that they 
correspond to two different \ values and for the system to undergo a transition 
from the chiral restored to the chiral broken phase it has to go through a high 
potential energy barrier. In a 3D plot of the energy per baryon number versus 
{p,x) the stable solutions correspond to two distinct "valleys", and the unstable 
solution mentioned at the end of the previous section corresponds to the top of the 
barrier between the two valleys. 



QUARK STARS 

In order to investigate the structure of stars we solved the Tolman-Oppenheimer- 
Volkoff (TOV) equation using the two EOS. Since EOS-I is identical to the one 
using a quadratic potential, it leads to stars that have the same phenomenology 
as the hybrid stars obtained by Drago et al. [15]: R ~ 10 — 12 km, a hadron crust 
and a mass M ~ 1 — 2M Q . At low densities, hadronization occurs and an hadronic 
equation of state should be used, replacing EOS-I. 

Since EOS-II saturates at a high density and, in addition, the system is not likely 
to undergo a transition to solution I, one should not perform any connection to the 
hadronic sector: the EOS-II alone generates a new family of strange quark stars. 
In Fig. 4 it is shown the mass-radius relation for different values of rj. These quark 
stars are smaller and denser in comparison with those resulting from EOS-I. For 
rj ~ 0.115 (and M = 1.7 GeV, yielding B 1 ^ ~ 0.195 GeV) one obtains a maximum 
radius R ~ 6 km and a corresponding mass M ~ 0.9M Q , which are the fitted radius 
and mass for the nearby compact object RX J 1856.5-3754 [20, 21]. According to our 
calculation, such star has a central density of 10p (Po is the nuclear matter density) 
and a central energy density e ~ 3xl0 15 g/cm 3 . At the edge, the density drops to 
5p and e ~ 1.35xl0 15 g/cm 3 . The ratio e/p remains approximately constant inside 
the star. The maximum period of the star, computed using the expression given 
in Ref. [24], is ~ 0.4 ms. From Fig. 4 one concludes that the mass-radius relation 
for these strange small stars mainly depends on the height of the local minimum 
of the x potential. 
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FIGURE 4. Mass versus radius for the pure quark stars (solution II) in the CDM model. The 
dot indicates the maximum radius star for r\ = 0.115. For the other model parameters see caption 
of Fig. 3. 

CONCLUSIONS 

Using a mean-field variational method we obtained two solutions for homogeneous 
strange matter in beta equilibrium using the CDM with a quartic potential, with 
its parameters fixed in the nucleon sector and to yield a reasonable bag constant. 

One solution is similar to the already known solution for quadratic potentials in 
the CDM, with massive quarks. In the other solution, quarks are almost massless. 
The EOS for both solutions are similar to those found recently in the framework 
of perturbative QCD. 

The pure quark stars emerging from the chiral symmetric solution are small 
and dense compact objects, showing a phenomenology compatible with the nearby 
X-ray source whose mass and radius was recently fitted [20, 21]. 
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